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ABSTRACT 



The T] — > 7r + 7r _ 7r°7 decay is discussed in the general context of Chiral Perturbation Theory 
(ChPT), assuming that the low-energy constants (counter-terms) are saturated by vector- 
meson resonances. The rj — ► 7r + 7r _ 7r 7 amplitude can be separated in two distinct pieces: 
the inner bremsstrahlung, A^ IB \ and the structure dependent (or direct emission), A^ SD \ 
amplitudes. The former - which essentially contains the same physics as A(rj — ► 7r + 7r~7r°) - is 
found to dominate over the second one - which looks more interesting from the ChPT point 
of view. 



Introduction 



At low energies, Chiral Perturbation Theory (ChPT) provides an accurate description 
of the strong and electroweak interactions of pseudoscalar mesons, P, in terms of a perturba- 
tive series expansion |]. At lowest order, ChPT is essentially equivalent to Current Algebra 
(CA) |2| but, at higher orders, loop effects appear restoring unitarity and giving rise to both 
finite corrections and divergences. The former are known to improve the lowest order (or 
CA) predictions, while the latter require the introduction of counterterms and destroy the 
renormalizability of the theory. The finite part of these counterterms, in turn, can be related 
to other successful aspects of hadron physics, such as the classical notion of Vector Meson 
Dominance (VMD). In this note we study the 77 — > 7r + 7r~7r°7 decay in the context of ChPT 
with VM dominated counterterms. 

The amplitude for the 77 — > 7r + 7r~7r°7 decay contains two distinct pieces receiving 
two distinct types of contributions, A(r) — > 7r + 7r~7T°7) = A^ IB ^ + A^ ^. The first piece - the 
inner bremsstrahlung amplitude A^ IB ^ - proceeds through the hadronic and G-parity violating 
i] — ► 7r + 7T~7r° decay, accompanied by the emission of an isovector (positive G-parity) photon. 
This part of the amplitude presents an infrared divergence and it is obviously dominated by 
soft photon emission. It can be deduced from Low's theorem and, therefore, A^ IB ^ essentially 
contains the same physics as the purely hadronic r] — ► 7r + 7r~7r° amplitude carefully discussed 
by Gasser and Leutwyler [|TJ . The second piece - the structure dependent (or direct emission) 
amplitude A^ SD ^> - looks a priori more interesting from the ChPT point of view. It involves an 
isoscalar (negative G-parity) photon and proceeds without the suppression due to G-parity 
violation. Therefore there is no obvious reason to assume that this G-parity conserving part 
of the amplitude, A^ SD \ has to be smaller than the first, G-parity violating one, 

A (IB) 

, except 

for the low-energy end of the photonic spectrum. 

The same conclusion was reached long ago by Singer || in his VMD calculation, as 
well as in the more detailed CA analysis of ref. (§], where A^ IB ^ was predicted to be larger 
than A^ SD ^ only for photon energies E y < 15 MeV, well below its maximum value E™ ax ~ 120 
MeV. These two and other old theoretical calculations for A^ SD ^ (based on CA and/or VMD) 
g |, § lead to 

. 11 = 0.23 %, 0.24 %, 0.28 %, (1) 

1 (77 — ► 7T U 77) 

although another similar (but oversimplified) analysis |J lead to a neglig ible A^^ and to 
smaller values for the ratio above. Combining the predictions ([I]) with present day data for 
the eta meson [0] one obtains 

T {SD \r] 7T + 7r"7r°7) - 2x 10~ 3 eV, B.RS sd \t] -> 7r + 7r"7r°7) ~ 1.6 x 10~ 6 . (2) 

2 



Up to now, this structure dependent part of the rj — > 7r + 7r 7r°7 decay has not been detected 
experimentally and only an upper limit - at 90 % of confidence level - is known H. 



r (5D) (r/ -> 7r + 7r-7r°7) < 6 x 10~ 4 T{r] -> a//) ~ 0.72 eV. 



(3) 



According to this discussion, the Saturne ^-factory and the Daphne 0-factory could in princi- 
ple produce enough etas to allow for detection and analysis of the structure dependent part of 
the 7] — > 7r + 7r~7r°7 amplitude. This prompted us to perform the corresponding ChPT calcu- 
lation following the lines of ref . || , where the analog (and similarly complicated) r\ — > 7r°77 
decay was considered. 

We will describe the nonet of pseudoscalar mesons P in terms of the SU (3) octet and 
singlet matrices 



\f2 V6 



71 



V 



7T 



K~ 



y/2 V6 



K° 



K° 

7E m J 



Pi 



(4) 



which appear in the ChPT lagrangian through the parametrization 

S = £ 8 £i = SxS 8 = exp (j (P 8 + Pi) J , 



(5) 



with / = 132 MeV |], [7|. Following the lines of ref. || and the reviews [T(|, we will adopt a 
rather simplified treatment of the singlet component of the rj in ChPT. More precisely, we will 
assume that the physical rj particle participates from nonet symmetry with an rj — v( mixing 
given by 



T) = cos 9 r/s — sin 9 rji 



[uu + dd — ss), sin #^—1/3, 



(6) 



as deduced from conventional rj — rj phenomenology and from more recent treatments in the 
ChPT context @. 

The lowest order lagrangian of ChPT (order two in particle four-momenta or masses, 
0{p 2 )) is 

f 



L 2 = J —tr{D^D^ + X £ f + X f S), 



(7) 



apart from an extra singlet mass term for 771 that should be added to account for the 
U(1)a problem, as discussed in [[TIJ]. The first term in (^) contains the covariant deriva- 
tive -D^S = <9 M X + ieAfj,[Q,Yl], with the photon field and the quark charge matrix Q 



[Q = diag(2/3, — 1/3, — 1/3)]. The non-derivative terms in Eq. (|7|), with x = = B A4, 
contain the quark mass matrix M. [M = diag(m n , ma, m s )] and lead to 

1 m 2 K m\ Am 2 K . . 



-B 



2 m u + m s m u + m d m d -m u 
where 

Am 2 K = (m 2 K0 - m 2 K+ ) QCD = (6.2 ± 0.5) x 10~ 3 GeV 2 (9) 

is an estimate for the QCD (or non-photonic) contribution to the squared K° — K + mass dif- 
ference. The latter plays a central role in both A(rj — > 7r + 7r _ 7r°) and the inner bremsstrahlung 
part of A(j] — > 7r + 7r _ 7r°7). Its numerical value in (^) is an average between the result [JTJ] 
Am 2 K = {m 2 K o ~ m \+ ~ m \® + m ^+) = 5-3 x 10 -3 GeV 2 , following from Dashen's theorem, and 
independent estimates |14]] (including improved versions of Dashen's theorem [|L5|| ) leading to 
Am 2 K in the range (6.5 - 7.0) x 10~ 3 GeV 2 . 

The next order lagrangian, 0(p 4 ), contains the Wess-Zumino term (the anomalous 
sector) and a series of ten (non-anomalous) counterterms identified and studied by Gasser 
and Leutwyler [I), 

L 4 = L wz + J2LiLf. (10) 

The only pieces of Lwz relevant for our purposes are the ones containing the anomalous 
PPPj and PPPPP couplings, i.e., 

e 



- x er#Ajr{Q%TldJ$dfiTtf - Qd^dJld^T) (11) 
2 



167T 

and 

rer^trtPdpPdvPdaPdpP), (12) 



157T 2 / 5 

respectively. The finite parts of the ten low-energy constants Lj, i = 1,...,10 are real and 
have been fixed by experimental data M . Alternatively, they can be approximatively deduced 



assuming that they are saturated by the exchange of known meson resonances [10|, [[J, HI 



thus justifying the phenomenological success of conventional VMD. Fixing the renormalization 
mass-scale around these resonance masses (/i = M p , for instance), the finite, renormalized 
values for Lj are small enough to justify the convergence of the perturbative series. This last 
remark obviously does not apply to L wz , generating anomalous processes with theoretically 
well-defined coupling strengths. 

The inner bremstrahlung amplitude A^ IB ^ 

As previously stated, this part of the amplitude for r](P) — ► n + (p + )n~ (p-)tt (po)j(q) 
can be related to the purely hadronic one, A(rj — > 7r + 7r~7r°), by using Low's theorem. Ei- 
ther through this theorem or by explicit calculation, one obtains the following, conveniently 



factorized expression 



A™ = - B(m f ~ mu) ± ( g± - ^] ( 1 + 2 m i~ 3P f + £/ + V + A (13) 



mt 



where the (infrared divergent) factor e (^j- — ) ( containing the photon momentum g and 
polarization e, comes from photon radiation by external pions, and the remaining factors 
correspond essentially to A(r) — > 7r + 7r~7r°) [|T|. 

The lowest order, 0(p 2 ), contribution proceeds through the tree level diagrams shown 
in Fig. la, b and c, with vertices of the L 2 lagrangian (|7|). At this lowest order, all r]—r]' mixing 
effects are ignored and the amplitude turns out to be given by eq. fll3D with U = V = W = 0. 
The two diagrams in Fig. la give the dominant contribution corresponding to the first term 
in the last parenthesis. The second term comes from the three diagrams in Fig. lb with a 
four-pseudoscalar vertex generated exclusively by the derivative part in L 2 . The remaining 
part of these three diagrams, generated by the massive couplings in L 2 , cancels with the 
contributions from Fig. lc. However, this lowest order - but otherwise exact - expression is 
expected to underestimate the r] — > 7r + 7r~7r°7 decay rate since the corresponding lowest order 
amplitude for A(j] — > 7r + 7r~7r°) is known to predict |jj a decay rate well below measurement 

0- 

The next order contribution to A^ IB \ 0(p 4 ), involves one-loop diagrams with vertices 
from (|7|) and tree-level counter-terms from the non-anomalous part J^LiL^ of (|10|). They 
lead to the U + V+W terms in eq ( |T3"D and their values will be taken from the detailed analysis 
of ref. [|TJ introducing two simplifying approximations. On the one hand, we will assume that 
the smallness of the available phase-space in both 7] — ► 7r + 7r~7r° and r] — ► 7r + 7r~7r°7 decays 
allows to approximate these corrections with the value at the center of the Dalitz plot for 
r] -> tt+tt-tt . With fi = 0.75 GeV ~ M p ~ Af w , this amounts to fix U + V = 0.39 - 0.03, 
coming from pion-loop effects, as discussed in detail in JTJ|. On the other hand, the dominant 
contribution from the finite part of the counter-terms is known to come Jl| from 7] — f]' 
mixing effects. In our nonet symmetry context with the conventional mixing angle, eq ([]), 
this amounts to take 1 + W ~ +y/2. Notice that this value corresponds to a value for W 
(where 7] — rj' mixing effects manifest) somewhat larger than that proposed in JL], but that 
our simplified treatment is essentially free from the criticisms recently rised by Leutwyler in 
ref. [ 18 1 . With these numerical values and eqs (|8|||) we obtain T(j] — > 7r + 7r~7r°) = 270 eV, 



in good agreement with experiment, Texp{v —* k + tt tt°) = 283 ± 27 eV [0. In spite of this 
agreement, we obviously do not claim that our simplified amplitude improves the original and 



more detailed ones in M and | lq| . We have simply achieved a successful parameterization for 



rj — > 7r + 7r 7T° from which we expect a reasonable prediction for A^ 1 ^ once inserted in eq (|13|). 



A similar treatment of the rj — ► 7r + 7r 7r° amplitude can be found in a recent analysis by Baur 



et al. 



For two different cuts in the photon energy, E 7 , we then obtain 

T iIB \r] -> 7T + 7r-7r°7) = 0.050 , £™ n = 50 MeV, 
r^^^vr+TT-A) = 0.76 eV , £f n = 10 MeV, (14) 

and the bremsstrahlung spectrum shown in Fig 2. Uncertainties affecting these predictions 
come mainly from the estimate in eq. @ and from neglected higher order corrections in 
ChPT, rather than from our simplified treatment of A(rj — > 7r + 7r~7r°). Globally, they should 
be expected to reach some 20 - 25 %. Our results ( |I~4"D are consistent with a previous analysis 
^0l once the two notations are unified. 



The Structure Dependent amplitude A^^ 

In contrast with the just discussed bremsstrahlung amplitude (0), which is propor- 
tional to the isospin violating factor — m u , ChPT predicts the existence of further, isospin 
conserving contributions to the global i] — > 7r + 7r~7r°7 amplitude. We now proceed to compute 
the dominant parts of these contributions and collect them into a structure dependent (or 
direct emission) amplitude A^ SD \ We will follow the ChPT analysis of 77 — > 7r°77 Q, whose 
amplitude A(r] —>■ 7r°77) - apart from the fact that it receives no inner bremsstrahlung con- 
tribution - is closely related to A( SD \ Indeed, in rj — > 7r°77 one photon is isoscalar (G = — ) 
and the other one is isovector (G — +). The former is the analog of the (isoscalar) photon in 
— > 7T + 7r 7T°7), while the latter - having the p° quantum numbers - plays the role of 
the 7r + 7r~ pair. Isospin symmetry allows to work in the good isospin limit, = m u , and to 
decompose A^ SD ^ in three terms by cyclically rotating the pion charge indeces 

A iSD)( v ^ ^-,,.07) = = A (SD \+, 0, -) + A (SD \0, -, +) + A (SD \-, +, 0). (15) 



The lowest order contribution to A^ SD ^ appears at order four and proceeds through the 
6 one-loop diagrams of Fig 3, containing two vertices from L 2 , eq (|7|). Part of this contribution 
is proportional to — m u and was already included in A^ IB \ but there is a second part - 
involving exclusively isospin conserving kaon loops - which belongs to A^ SD \ This part has to 
be finite since no suitable counterterms are available in ( JTOj ) for uid = m u . One finds indeed 
the finite result 

A[ 4 7 ) (+,0,-) = 6 — (6Qq + 3Q 2 - 4m 2 K )[(ep + )(qp-) - (ep-)(qp + )]I(m 2 K , Q\ Qq) (16) 
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with Q = p + + p_ and 

-1 rl-X 



I(m 2 K ,Q 2 ,Qq) = I dx I dy— — — — (17) 

Jo Jo m 2 K - 2Qqxy - Q 2 y{l - y) 

By itself, this order four contribution leads to 

rjf) D) (r/ -> 7T + 7r-7r°7) = 0.88 x 10~ 7 eV, (18) 

well below the old estimates ([!]) and (@), and suggesting that this is not the dominant contri- 
bution to A^^. The same situation was found when analizing 0(p 4 ) kaon loops for rj — > 7r°77 
in ChPT [[J . Similarly, the smallness of this lowest order contribution seems to be a reminis- 
cence of the vanishing of the old CA estimate of ref . || . 

At next order, 0(p 6 ), one has contributions coming from tree level (or counterterms) , 
from one loop and from two loops. The former belong to L$ and their finite part will be 
obtained from saturation with vector mesons. For further reference, we compute the full 
VMD amplitude for rj — > 7r + 7r~7T°7, which proceeds through the two diagrams shown in Fig 
4 (apart from rotations of pion indexes). Ignoring negligible effects from T pu finite widths, 
one obtains 



ASD) , 



4vr 4 / 2 (M2 - Q 2 ) 



P P' 



P p P w i(ep + ) ((Pp-)(qp ) - (Pp )(qp-))] + 

(19) 



(ep+) {{Pp-){gpo) - {Ppo){gp-)) - ^{qpo) {{tp+){qp-) - (ep_)(gp+)) 



with 



PP = 1 ( 1 + 1 + 1 



" ~ " M 2 -(P-q) 2K M 2 -Q 2 M 2 p -Q% M 2 -Q 2 _ 



P n P' 



11 11 11 



M 2 - Q 2 M 2 - (P - Q) 2 M 2 -Q 2 + M 2 -(P-Q + ) 2 M 2 - Q 2 _ M 2 - (P - Q_) 

and Q± = p± + p . The coupling constants are such that M 2 ~ M 2 ~ 2f 2 g 2 and g = 4.2, 
as discussed in refs.|9|, [12], [i~3| . The part of this VMD contribution which corresponds to 
the Lq counter-terms in ChPT can be simply obtained from ([T^) by expanding the VM 
propagators 1/(M V — K 2 ) = l/My + K 2 /My... and retaining only the first term. This leads 
to FvMD(G){ r l ~^ 7r + 7r _ 7r°7) = 0.42 x 10 -4 eV, which turns out to be much larger than the 
lower order estimate ([18]). Since in ChPT loop corrections at order six are expected to be only 
a fraction of the corresponding loop corrections at order four, eq (|TJ), we can safely conclude 
that the full ChPT prediction at order six is drastically dominated by the VM saturated 
counter-terms, i.e., 

r% D) (v - vr+TT- A) ^ rv D MD(6)(V - tt+tt-A) = 0.42 x 10- 4 eV (20) 
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The same dominance of 0(p 6 ) counterterms over the corresponding loops was also observed 
in the case of A(r) — > 7r°77) ||. 

At order 0(p 8 ) more counterterms appear and a new type of loop-correction becomes 
potentially important. The contributions of resonance dominated counterterms can simply 
be obtained by expanding the full VMD amplitude (|l9|) and retaining only 0{p 8 ) terms. By 
itself, this leads to a non negligible order p 8 correction, namely, ^vmd(8)(v ~ > 7r + 7r~7r°7) = 
0.094 x 10~ 4 eV representing a non negligible increase to Eq . (|20"D . The new type of loop effects 
looks a priori more interesting. Taking two vertices from the anomalous Lwz (one from flllD 
and another from (|12"D) one obtains a non-anomalous one-loop correction of order p 8 , which 
does not vanish only for kaon loops. This "doubly-anomalous" loop contributions leads to a 
decay rate of 4.2 x 10~ 7 eV, i.e., of the same order as Tu^\tj — > 7r + 7r~7r°7) and well bellow 
the order eight counterterm. Again, we find that kaon loops for A^ SD \ri — > 7r + 7r _ 7r°7) closely 
follows the same pattern as kaon loops for A(rj — > 7r°77). The whole order eight contribution 
alone is therefore also dominated by counterterms, 

r[ 8 5 )% - A) -r^ (8) (?7^7r + 7r-A) = 0.094 x 10" 4 eV (21) 

All this implies that in a ChPT context with resonance saturated counterterms, the 
whole structure dependent amplitude for 77 — > 7r + 7r~7r°7, A( SD \ is strongly dominated by 
the contributions of these counterterms and should essentially be given by the full VMD 
amplitude (0). Such an "all-order" estimate leads to 

T SD (r] -> 7r + 7r"7T°7) ~ T V md(v -> vr + 7r"7r°7) = 1.4 x 10 -4 eV (22) 

The corresponding photonic spectrum is shown in Fig. 2. 

This structure dependent part of the amplitude is the one that should be compared 
to older estimates, although none of those treatments coincides precisely with ours. From 
eq (|22|) and the result Tvmd(v — » 7r°77) = 0.31 eV of ref. 0, one obtains Tvmd(v — > 
7r + 7r~7r°7)/TyA/£)(?7 — > 7r°77) ~ 0.44 x 10~ 3 somewhat below the old estimates (|1|) and above 
the vanishing prediction of ref. 0. The shape of our photonic spectrum coincides with the 
old prediction by Singer || who used a simplified version of VMD not far from ours. 

Our final ChPT predictions for the whole rj — > 7r + 7r~7r°7 amplitude can be obtained 
from the sum of A^ IB \ eq flilf), and the full VMD amplitude, eq (fl9|) . The contribution of 
the latter - in modulus plus interference with A^ IB ^ - has negligible effects for soft photons 
and represents a minor increase (somewhat below 1 %) of the decay rate in the higher half of 
the photonic spectrum as shown in Fig. 2. The integrated width remains thus unaffected as 
in eq (0). 

8 



Conclusions 



From the preceding and somewhat intricate analysis of the rj — > 7r + 7r~7r°7 decay, a 
definite two-fold conclusion emerges, namely, that the whole amplitude is strongly dominated 
by inner bremsstrahlung and that it is expected to be large enough to allow for detection in 
a near future. This is a useful conclusion in that it contradicts (and presumably corrects) 
older results (see refs. and it will furnish a clear-cut test for (and presumably confirm) 

ChPT. But it is also a deceptive conclusion reducing most of the 7] — > 7r + 7r~7r°7 dynamics to 
that of 7] — > 7T + 7r~7r°, already studied both theoretically |IJ and experimentally ]7] (although 
not free of uncertainties [|], [18], ^1[). Only an extremely sensitive and dedicated experiment, 



and an improved theoretical treatment of A^ 13 " 1 could allow to extract the structure dependent 
part of the amplitude which contains the genuinely new effects of ChPT for rj — > 7r + 7r~7r°7. 
In this case, the values of counterterms and the hypothesis of their resonance saturation - 
rather than chiral loop effects - will be tested. Stated otherwise, the experimental analysis 
of 7] — > 7r + 7r~7r°7 can represent an excellent confirmation of ChPT if the predicted spectrum 
(largely dominated by bremsstrahlung) is observed, but it can hardly be useful to improve 
our knowledge on other aspects of this theory . 
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b) 



c) 




Figure 1: Diagrams contributing to A IB (r] — ► 7r + 7r 7t°7) at lowest order, 0(p 2 ), in ChPT. 
The cross in each diagram indicates a coupling proportional to — m u . 
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Figure 2: Spectrum of the photon energies, E 1) as predicted by ChPT for the 77 — > 7r + 7r~7r°7 
decay. The solid line corresponds to the inner bremsstrahlung contribution. The dashed line 
is the VMD contribution alone. The dotted line is their interference. 
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